A note on the geometric phase in adiabatic approximation 
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Abstract 

The adiabatic theorem shows that the instantaneous eigenstate is a good approximation of the 
exact solution for a quantum system in adiabatic evolution. One may therefore expect that the 
geometric phase calculated by using the eigenstate should be also a good approximation of exact 
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\q ' geometric phase. However, we find that the former phase may differ appreciably from the latter if 
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tj" ■ the evolution time is large enough. 
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The existence of Berry phase was demonstrated in 1983 [jj. In essence, Berry showed 
that a quantal system in an instantaneous eigenstate of its Hamiltonian, when slowly trans- 
ported round a circuit by varying parameters in its Hamiltonian, will acquire a geometric 
phase factor in addition to the familiar dynamical phase factor. Here he made explicit use 
of the adiabatic theorem which states that if a quantum system with a time- dependent 
nondegenerate Hamiltonian H(t) is initially in n-th instantaneous eigenstate of H(0), and 
if Hit) evolves slowly enough, then the state of the system at time t will remain in the n-th 
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instantaneous eigenstate of H(i) up to a multiplicative phase factor [2, |3j, |4| . This discovery 
has prompted a myriad of activities that has led to its generalization. Indeed, Anandan 
and AharonovJ5| extended Berry's result to nonadiabatic cyclic evolutions. J. Samuel and 
R. Bhandari[6] further generalized the geometric phase to noncyclic and nonunitary evo- 
lution. Essentially, the original definition of geometric phase found by Berry for quantum 
systems under conditions of adiabatic, unitary and cyclic evolutions has been generalized to 
nonadiabatic, nonunitary and noncyclic evolutions [21, |SJ. Other generalizations include geo- 
metric phases for mixed states, geometric phase for open systems and off-diagonal geometric 
phases^ bj| U, |jj, U , Ili|l5 1 

Recently, some authors [ly, |l7|, |l8( noted that one must be careful in employing the adi- 
abatic theorem. Indeed, Marzlin and Sanders demonstrated that a perfunctory application 
of the theorem may lead to an inconsistency |l6J; and Pati and Rajagopal observed that 
the Berry phase might be zero if some terms like (-Em|-£/fc)|lS| are incorrectly dropped and 
they pointed out that these terms contribute to the Berry phase. These works led us to ask 
whether there is any significant difference between the approximate phase evaluated under 
the assumption of adiabacity and the exact geometric phase even if the adiabatic theorem 
is carefully used. 

In this note, we compare the geometric phase calculated under the assumption of adiabatic 
theorem with the exact geometric phase. In examining a well-studied example of a spin-half 
system in the presence of a rotating magnetic field, we find that the phase calculated by 
using the adiabatic approximation may differ appreciably from the exact geometric phase if 
the evolution time is large enough. We further show that the difference is inherently present 
in general systems. 

Consider a closed iV-dimensional quantum system S with the Hamiltonian H(t). The 



instantaneous nondegenerate eigenvalues and orthonormal eigenstates of H(t), denoted as 
E m (t) and \E m {t)) respectively, are defined by 

H(t)\E n (t)) = E m (t)\E m (t)), m = l,...,iV, (1) 

Suppose the system is initially in one of the eigenstates |^(0)) = 1^(0)) and its time 
dependent state is denoted by \ip(t)). The geometric phase 7(7") obtained by the system 
after experiencing an evolution from t = to t = t can be calculated by using the well- 
known formula [3], 

T (r) = arg^(0)|^(r)> +1 [ {^{t)\^{t))dt, (2) 

Jo 

which is valid for both nonadiabatic as well as adiabatic evolutions. 

Specially, if the Hamiltonian H(t) changes slowly and E m (t) obey the adiabatic con- 
straints 

(E m (t)\E k (t)) 



E m (t)-E k (t) « 1 > m ^ k > ( 3 ) 

according to the adiabatic theorem, the state of the system at time t will be 

m))*j a(t) \E k (t)), (4) 

where a{t) is a phase parameter. Then, by substituting Eq. (JIJ) into Eq. (J2J), we get the 
geometric phase in the using of adiabatic theorem as 

T (r) = arg(£ fc (0)|£ fc (r)) + i f ' (E k (t)\E k (t))dt, (5) 

■J 

which may be regarded as the generalized noncyclic Berry phase, extended from Berry's 



original definition of cyclic Berry phase 
the form of 



m 



It is worth noting that the Berry phase takes 



l{r)=i (E k {t)\E k {t))dt (6) 
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if and only if a.Tg(E k (0)\E k (r)) = 0. We would like to stress that Eq. (jSJ) is valid 

for both cyclic and noncyclic adiabatic evolutions, while Eq. (JHJ) is valid only when 

arg(£ fc (0)|£ fc (r)> = 0. 

It is instructive to ask whether there is any appreciable difference between the approxi- 
mate geometric phase and the exact geometric phase for the same system under adiabatic 



evolution. In other words, we would like to compare the phase calculated using Eq. (JSJ), 
where the adiabatic theorem is applied, with the phase defined by Eq. (J2J with the exact 
state. The widely accepted notion is that there is little or no difference between the two. As- 
tonishingly, our result shows that the difference is appreciable if the observing time interval 
is large enough. 



Let us first consider a well- 
which was also considered in 



mown model; a spin-half particle in a rotating magnetic field, 
The Hamiltonian of the model can be written simply as 



H(t) = (a x sin 9 cos uot + a y sin 9 sin ut + o z cos ( 

uj I cos 9 sin9e~ luJt 
2 ^sinfle*"' -cos 9 

where ujq is a time-independent parameter defined by the magnetic moment of the spin and 
the intensity of external magnetic field, 10 is the rotating frequency of the magnetic field and 
<jj, i = x,y,z, are Pauli matrices. The instantaneous eigenvalues and eigenstates of H{t) 
are 

E i = y> E 2 = — g-, (8) 

and 



\Ei(t)) 
\E 2 (t)) 



e -iut/2 gin | 

_ e iu J t/2 cog ^ 



e iuJt/2 sin 2 




(9) 



respectively. The adiabatic constraints Q are satisfied as long as 

u ^>ujsm9. (10) 

Suppose that the system is initially in the state \Ei(0)). At time t, according to the adia- 
batic theorem, it will be in the instantaneous eigenstate \Ei(t)) up to a phase factor. The 
geometric phase for the evolution \Ei(t)), t G [0,r], can be calculated by Eq. (J3J), and we 
have 

j(t) = arg I cos — + i cos 9 sin — ) — -lot cos 9. (11) 



When r = T = 2tc/uj, the cyclic geometric phase reduces to 7(T) = tt(1 — cos^). 

Next, we evaluate the geometric phase of the above model without any approximation. 
Here, we solve the Schrodinger equation 

/ 
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UJ Q 



if t m)) - -y 



cos 9 sin i 



-iuit 



sin i 



Auit 



cos 9 



\m), 



with the initial condition \ip(0)) = |£'i(0)), and obtain the exact solution, 

m)) = a(t)\Ei(t))+b(t)\E 2 (t)), 



where \Ei(t)), \E2(t)) are given by Eq. Q and a(t), b(t) are given as 

a{t) 

m 



Ujt UJq + UJ cos 9 uot 
cos i sin — , 

2 uj 2 ' 



uj sin 9 iot 
% sin — , 

u 2 ' 



(12) 



(13) 



(14) 



with uj = Jul + uj 2 + 2uj uj cos 9. Substituting Eq. (|T3*j) into Eq. (J2J), we obtain the geomet- 
ric phase, 



7(t) 



arg 



UJT UJT UJ + UJq COS 9 UJT UJT 

cos — cos 1 sin — sin — 

2 2 W 2 2 



( UJT UJT UJQ + UJCOs9 UJT UJT^ 

+ i cos 9 cos — sin sin — cos — 

\ 2 2 uj 2 2 / 

ujquj 2 sin 2 9 __ uj q t uj uj 2 sin 2 9 
H :r^3 SID. UT H — ; ;r-s 7" 



(15) 



2cJ 3 2 2cJ 2 

Eq. ()15|) gives the exact geometric phase of the system without any approximation. It is 
valid both for adiabatic evolution as well as for nonadiabatic evolution. Now, if we invoke 
the adiabatic condition, uj q ^> uj sin 9, we have 



/ \ UJT . . UJT 

7(r) ~ arg I cos h*cos0sin — 

1 ujquj 2 sin 2 9 

— -UJT COS 9 7, T. 

2 2cJ 2 



(16) 



If (uj uj 2 t sin 2 9) /(2uj 2 ) <^C 1, we further have 



I \ ( UJT . . UJT 

7(t) ~ arg I cos — + z cos & sin — 

— -ujt cos 9, 



(17) 



which is the same with Eq. (Illj) . From the above, we find that both cyclic phase 7(27? /uj) 
or noncyclic phase 7(7") do correctly describe the geometric phase in adiabatic evolution if 
the time r is not too large, say r e (0, 2tt /uj). We also note that the phase defined by Eq. 
(JHJ) may not give the correct geometric phase even in adiabatic evolution if the time r is 
large enough, in which case, the difference between the exact and the approximate value, 
£7 = — (ujouj 2 t sm 2 9)/(2uJ 2 ), can become significant. To elucidate this, we note that with 
uj 2 ~ Uq + 2ujquj cos 6, Sj can be written as 

uj 2 sin 2 6 



(hj 



2{uq + 2uj cos 9) 

(18) 



where s = to> 2 sin 2 9/2{ujq + 2ujcos9). Under the adiabatic condition, s is a small quantity. 
So, £7 is a small quantity in general. However, £7 may become an appreciable quantity if 
t is large enough. Although in real experiments r is usually taken to be one or two cyclic 
periods, in theory it may be much larger than one. Note that the adiabatic theorem is 
valid even for r 3> 2tt/uj as long as the adiabatic constraints (jl(J|) is satisfied. One may 
verify that \ip(t)} in Eq. (fT3|) always reduces to \Ei(t)) in Eq. Q up to a phase factor, 
independently of time, as long as the condition (fTTH) is satisfied. More accurately, one may 
compare the rays in the projeted Hilbert space. Let Ap = p^, — p E , where p^ = |\l/(t))(\I/(t)| 
and pE = \Eiit)){Ei{t)\. One obtains 

Ap= -(^) 2 sin 2 f (\E l (t))(E 1 (t)\ - \E 2 (t))(E 2 (t)\) 
u 2 

, (cun + uj cos 6)u) sin 6 . , uJt .ujsm.6 cut . ZJt,.^,,.,^,,. 

sin 2 — + 1 _ cos — sin— )\E 1 (t)){E 2 (t)\ 



uJ 2 "2 Lj 2 2 

(uj + uj cos 9)uj sin 9 . 9 UJt .ujsm9 UJt . Tut. .,_,....,_,. .. ._. 
sin 2 z _ cos — sin— ) \E 2 {t)) {E x (t) \, (19) 



v uj 2 2 uj ' 2 ~ ' 2 

and will find Ap — > 0, in the condition where uj$ ^> uosmO. We may therefore observe 
that the eigenstate \Ei{t)) is always a good approximation of the exact solution \ip(t)) for 
arbitrary times t, but the geometric phase calculated by using the eigenstate may not a good 
approximation of the geometric phase calculated by using the exact solution if t 3> 2k /uj. 

We now extend the above analysis to a general iV-dimensional quantum system with 
Hamiltonian H(t). If the system is initially in the state \ip{0)) = 1^(0)) then its exact state 
at time t can be written as 

m)) = Y. C m(t)\E m (t)), (20) 



where C m (t) satisfy Xi m Cm(i)C m (t) = 1, C m (0) = 5 mk . Without any loss, we may further 
let 

C m (t) = e m (t) + 6 mk e ta{t \ (21) 

where the parameters e m {t) satisfy e m (0) = and 

J2e* m (t)e m (t) = s k (t)e-^ - effle**®, (22) 

m 

and a(t) is defined by Eq. (jlj. Substituting Eqs. (J2UJ) and (J2TJ) into Eq. (J2J), we obtain the 
geometric phase as (see Appendix for details) 

7 (r) = arg(E k {0)\E k (r))+i f {E k {t)\E k {t))dt 
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+A 7 , (23) 



where 



1 + V£ () (E k (0)\Em(T)) (T) 

4- ml V(o)i^(r)> e 



A7 = arg 

+iU k (T)e-^^ + J2J o £* m de m ) 

-2Re ( T e k {t)a{t)e~ ia{t) dt\ 

-2Im( £j2 £ rn(t)(E k (t)\E m (t))e- ia ®dt) 

- Im [f Q T. e Ut)en{t)(E m {t)\E n {t))dt\ . (24) 



Eq. ()23|) with Eq. (J24J) represents the exact geometric phase without any approximation. 
Notice that the second term on the right side of Eq. (J24}1 is also real, as a consequence of 
Eq. (EU). 

Comparing the phase defined by Eq. (jHJ) with the exact geometric phase, we find that 
A7 is just the difference between them. For the quantum system in the adiabatic evolution, 
\ip(t)) ~ e ia ^\E k (t)), we must have |e m (£)| ^ 1- A7 approximates to zero if the time r is 
not too large, and in this case, the approximate phase does properly describe the geometric 
phase. However, if the time r is large enough, such that A7 is not a small quantity, the 
phase defined by Eq. (J3J may be very different from the exact geometric phase. This is 
because the last three terms on the right of Eq. (|2*4*|) may become appreciably large as 



the time progresses unlike the first two terms which can be neglected under the adiabatic 
approximation. It is important to notice that the property of phase is very different from 
that of other physical observables. For example, if the number (2 10 + 1)tt is the value of 
energy, then it is appropriate to take 2 l0 ir as its approximation; but we cannot do this if it 
is the value of phase. In other words, an infinitesimal error over a small time interval may 
accumulate to a measurable quantity as the time becomes large. 

In conclusion, we have examined the difference between the geometric phase calculated in 
adiabatic approximation and its exact counterpart. Here we have showed that the difference 
can become appreciable over long time intervals. This result, in our opinion, is significant, as 
it debunks the commonly held notion that the two phases are no appreciable difference under 
the adiabatic approximation. It alerts us that although the adiabatic approximation of wave 
function is valid, its further application needs to be earful especially in the calculation of 
geometric phase. This theoretical analysis, however, may not affect existing experimental 
results since the evolution time in the existing experiments is usually short; within one or 
two cyclic periods. We also would like to stress that i J T (Ek(t)\Ek(t)) is not generally gauge 
invariant. Indeed, even in the cyclic adiabatic case, the presence of arbitrary phases in the 
eigenstates of a time-dependent Hamiltonian may introduce non-trivial contribution to the 
above term. It should be noted that if the eigenstate satisfies arg(£ , fc(0)|-E'fc(r)) = 0, then Eq. 
(JBJ) gives the Berry phase. In general, however, for a system undergoing adiabatic evolution 
one should use Eq. (J5J) to calculate the Berry phase. 

We thank A. K. Pati for useful remarks. The work was supported by NUS Research 
Grant No. R-144-000-071-305. 

Appendix 

In the appendix, we outline some of steps involved in obtaining the expressions (|23|) and 

(123). 

Substituting Eqs. (j2"Uj) and (}2"T]) into arg(^(0)|^(r)), by directly calculating we have 
arg(V(0)|V(r)> 



arg 
arg 
arg 
arg 



(E k (0)\J2C m (r)\E n (r)} 

m 

E { £ m(r) + 8 mk e ia ®) (E k (0)\E m (r)) 

m 

e ta ^(E k (0)\E k (r))+Y,^m(r)(E k (0)\E m (r)) 



e ia ^(E k (0)\E k (r))(l + J2e m (r) 



(E k (0)\E m (r)) 

(E k (0)\E k (r))' 



-ia(r) 



a(r) + arg(£ fc (0)|£ fc (r))+arg 



1 + V£ , (E k (0)\E m (r)) ia{T) 
^ Z } (E k (0)\E k (r)) e 



(25) 



where the condition (E k (0)\E k (r)) ^ has been used [20]. 

Substituting Eqs. (J20I) and (J2IJ) into i Jq (ip(t)\ip(t))dt, we have 



i / m)m))dt 







d 



* l Q E C *m(t)(E m (t)\ Jt r£ C n(t)\En(t))\ dt 



^C(^n({)(^(i)|K(<)> + EC(i)^(t){^WI^(i)> 

m,n m,n 

E (C(*) + ^e- ia (t) ) (e m (t) + i6 mk a(t)e ia ® 



dt 



dt 



+ E ^* m (t) + 5 mk e- ia ®) (e n (t) + 5 nk e ia ®) (E m (t)\E n (t)) 

it) + (E k (t)\E k (t)) + ie* k (t)a(t)e ia ® + e k {t)e' ia ^ + 1 £e* m (t)e m (t) 

m 

+ EC(*)e iQ(t) (^(t)|4(t)}H-E^(t)e-^( J B fc (t)|K(t)) 



za 



+ E4(iKW(^(t)IKW) 



di 



(26) 



Noting that 

/" T v^ v^ f £m{ - T) 

i E £ m( t )^( t ) rft = *E / £ m^ 

jo ™ tf Jo 



f£m(r) 



z r [i4(t)d(t)e iQ(4) + e k {t)e- ia{t) ] dt = ie k (r)e- ia(T) - 2Re ( f" e k {t)a{t)e-' ia{t) dt\ 
i£j2 £ *m( t ) e n(t)(E m (t)\E n (t))dt = -Im ( [ Y d e* m {t)e n {t){E m {t)\E n {t))dt ) , 



m,n 



and 



E £ m(*)e fa(t) <S m (t)|4(*)) +Y, £ n(t)e~ ia V{E k (t)\E n (t)) ) dt 



= -21m [ T Y,e m (t)(E k (t)\E m (t))e- ia ^dt) , (27) 

we obtain 

i [ T (m\m)dt 



= -a(r)+tj o (E k (t)\E k (t)}dt + t(e k (r)e- ia ^+J2j o e* m de m \ 

- 2Re ^£ e k (t)a(t)e- ia ^d?j -2lm(j^^rn{t){E k {t)\E m {t))e~ ia ^dt\ 
-Im ( f T Y d e* m {t)e n {t){E m {t)\E n {t))dt\ . (28) 

\^° ran / 

Taking Eqs. (j2*5Jl and (|2*8j) into (0), the expressions (|2*Hj) and ()24j) are immediately obtained. 
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